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We study the surface critical behavior of semi-infinite systems be- 
longing to the bulk universality class of the Ising model. Special 
attention is paid to the local behavior of experimentally relevant 
quantities such as the order parameter and the correlation func- 
tion in the crossover regimes between different surface universality 
classes, where the surface field hi and the temperature enhance- 
ment c can induce additional macroscopic length scales. Starting 
from the field-theoretical (j)'^ model and employing renormalization- 
group improved perturbation theory (e expansion), explicit results 
for the local behavior of the correlation and structure function (0- 
loop) and the order parameter (1-loop) are derived. Supplementing 
earlier studies that focussed on the special transition, we here pay 
particular attention to the situation where a large c suppresses 
the tendency to order in the surface (ordinary transition) but a 
surface field hi generates a small surface magnetization mi. Our 
results are in good agreement with recent phenomenological con- 
siderations and Monte Carlo studies devoted to similar questions 
and, combined with the latter, provide a much more detailed un- 
derstanding of the local properties of systems with weak surface 
fields. 

Key words: Field theory, renormalization group, surface critical 
phenomena, critical adsorption 



Preprint submitted to Elsevier Preprint 



1 February 2008 



1 Introduction 



Surface critical phenomena are subject of increased current experimental in- 
terest. After the confirmation of some of the theoretical predictions [1-4] in 
X-ray scattering experiments with FeAl [5], the more recent efforts focussed 
on binary mixtures near their consolute point [6-8] and near-critical fluids [9]. 
While for instance the diffuse scattering of (evanescent) X-rays is governed 
by the two-point correlation function, in the experiments on fluids the order 
parameter (OP), the concentration difference in fluid mixtures or the density 
difference between liquid and gaseous phase in single-component fluids, plays 
a central role. For instance the reflectivity and ellipticity measured in light- 
scattering experiments are directly related to the OP profile[10,ll]. Hence, a 
precise quantitative information about the local behavior of these quantities 
is required to interpret the experimental data. 

In the framework of continuum field theory such as the 0^ model (belonging 
to the universahty class as the Ising model) the surface influence is taken into 
account by additional fields hke the surface magnetic field hi and the local 
temperature perturbation c at 2; = 0. The latter can be related to the surface 
enhancement of the spin-spin coupling in lattice models [2]. At the bulk critical 
temperature, r = {T — T^)/T^ = 0, the tendency to order near the surface can 
be reduced (c > 0), increased (c < 0), or, as a third possibility, the surface 
can be critical as well (corresponding to a particular value of c). As a result, 
each bulk universality class in general divides into several distinct surface 
universality classes, corresponding to the fixed points of the renormalization- 
group flow and labelled as ordinary (c — > 00), extraordinary (c — > —00), and 
special transition (c — c*). 

While a very well-developed theory for the individual universality classes ex- 
ists, the picture in the crossover regions between the fixed points is less com- 
plete. In particular in view of some of the experiments, it may not be justified 
to consider the system to be at a fixed point, and so a detailed understand- 
ing of the crossover regions is of special importance. Additionally the physics 
away from the fixed points is much richer and more interesting as certain 
length scales emerge, may become macroscopic, and compete with the bulk 
correlation length, whereas with the surface fields at their fixed-point values 
these scales are typically zero or infinity. 

Near the special transition these phenomena were studied already some time 
ago[12,13]. It was realized that hi gives rise to a length scale Igp ^ hi'^ ^ 
and the singular behavior of thermodynamic quantities goes through different 
universal regimes. This is best illustrated by the spatial dependence of the OP 
m{z). For distances large compared to all length scales it was known that m{z) 
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should decay as ~ z'^^" [1], where the exponent P/u is the scaling dimension 
of the field (f). This power law describes the asymptotic long-distance decay at 
when the symmetry in the surface is broken spontaneously or explicitly, i.e. 
at the extraordinary or normal transition, respectively[14,15]. For instance for 
the d — 3 Ising model ~ 0.52 [16]. For z <^ Igp (but still much larger than 
microscopic scales), on the other hand, the result found in Rcfs. [12,13] was 
that ■m(z) should behave as ~ z^^^ -z^)/'^. Again, in the three-dimensional Ising 
model the Monte Carlo literature value for this combination of exponents is 
{Pi" — ~ —0.15 [17], whereas in mean- field (MF) theory it is zero. Thus, 
at the special transition, if fluctuation are taken into account, the decay of 
■m{z) is still monotonic, however, it is described by different power laws for 
2; — > and z — > 00. 

A similar phenomenon was recently pointed out for the case of large c, i.e. 
a situation near the ordinary transition [18]. For c » c* the scahng held is 
given by hi = hi/c^ [2] (where y is a positive exponent to be discussed below 

in more detail) and the length scale induced by this fleld is lord ~ hi ^ . 
The short- distance behavior of the magnetization is given by m ~ with 
K = (A°'"'^ — f3)/v. This time the value of the short-distance exponent, /t, is in 
general positive, such that the OP turns out to be a non-monotonic function 
of z. In the Ising model k ~ 0.21, whereas in MF theory, hke at the special 
transition, one has k — At z ^ lord again the crossover to the "normal" 
behavior z~^^'^ takes place. The near-surface growth of order at the ordinary 
transition was recently also corroborated by Monte Carlo simulations [19]. 

Moreover, it was pointed out that in a system with c ^ c* and non-vanishing 
hi , it might be the generic case that length scale Zord becomes macroscopic and 
comparable or even larger than the bulk correlation length [20] . The behavior 
of thermodynamic quantities near the surface is then governed by exponents of 
the ordinary transition. As a consequence, for ^ ~ lord the amount of adsorbed 
order as a function of r is described by the power law t^'^T [20], in very 
good agreement with experiments [9,8]. Only much closer to the critical point 
a crossover to t^^'^ takes place, as expected for the normal transition and 
assumed in previous studies. 

The main purpose of our present work is a careful derivation of the OP scal- 
ing function for the crossover between ordinary and normal transition in the 
framework of renormalization-group improved perturbation theory. While the 
respective studies in Ref. [18] focussed on the near-surface behavior, and the 
Monte Carlo simulations of Ref. [19] were hampered by strong finite-size effects 
and did not yield OP profiles suitable for the comparison with experimental 
data, the results presented below give to one-loop order the OP profiles for the 
semi-infinite system. Further we present MF results for the structure function. 
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relevant for scattering of X-rays for example, in the regime between ordinary 
and normal transition. Eventually we discuss the crossover between special 
and ordinary transition that (while probably less important in view of exper- 
iments) is an interesting and technically demanding problem in perturbation 
theory. 

The rest of this paper is organized as follows: In Sec. 2 we discuss the behavior 
of thermodynamic quantities from the viewpoint of a scaling analysis, give 
an heuristic argument for the growth of order at the ordinary transition, and 
compare with the situation in critical dynamics where analogous phenomena 
occur after temperature quenches. In Sec. 3, as the building blocks for the 
pcrturbative calculation, we present the results of the MF (zero-loop) theory. 
Especially the correlation function for general c and hi is discussed in some 
detail, and the structure function for critical diffuse scattering is calculated 
numerically from that result. In Sec. 4 the OP profiles are calculated to one- 
loop order, and special emphasis is put on the crossovers between ordinary 
and normal transition and special and ordinary transition, respectively. 



2 Background 

2. 1 Scaling Analysis for the Order Parameter 

In this section we study the behavior of thermodynamic variables with the help 
of a phenomenological scaling analysis. As the most instructive example, we 
discuss the local behavior OP for different surface universality classes and the 
various crossovers. Certain aspects of the scaling behavior of the correlation 
function will be treated in Sec. 3. 

Since the parameter c has to be fine tuned in addition to r, from the viewpoint 
of the experimentalist the special transition is certainly the most exotic among 
the surface universahty classes. However, from the viewpoint of the scahng 
analysis it is the most straighforward case. Near the special transition hi and 
c = c — c* are the linear scaling fields pertaining to the surface. In the critical 
regime thermodynamic quantities are described by homogeneous functions of 
the scaling fields. Let us consider the OP for small hi and c > 0. Its behavior 
under rescaling of distances should be described by 

m{z, T, /ii, c) ~ 6"^/'^ m{z 6" \ Tb"'" , hi c b^'") , (2.1) 

where all exponents have their standard meaning[l]. In general the surface 
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exponent Ai occurring in (2.1) has different values for different surface uni- 
versality classes[l], and therefore it has been marked by 'sp' (for special). 

Removing the arbitrary rescaling parameter b in Eq. (2.1) by setting it ~ z, 
one obtains the scaling form of the magnetization 

m{z, T, h) - z-^/''M{z/^, z/hp, z/Q , (2.2) 

where 

hp ^ hi"^^'" and l^r^c""/^ (2.3) 

are the length scales determined by hi and c, respectively. The other length 
scale pertinent to the semi-infinite system and occurring in (2.2) is the bulk 
correlation length ^ ~ t~'^. 

In order to discuss the various asymptotic cases, let us set r = for simplicity. 
In other words, we assume for the moment that the bulk correlation length is 
much larger than any length scale set by the surface fields, in which limit the 
scahng function A4 in (2.2) becomes 

M{z/^, Z/Isp, Z/Q ^ Mcrit{z/lsp, z/lc) ■ (2.4) 

Further below we will discuss modifications due to finite ^. Let us first consider 
the case c = where both bulk and surface are in the critical state at T^. Then 
the only remaining length scale is Igp and the OP profile can be written in the 
scaling form 

m{z, hi) ~ z-f^/" Msp{z/lsp) . (2.5) 

As mentioned in the Introduction, for z — > oo the magnetization decays as 
~ z~^^^ and, thus, M.sp{0 should approach a constant for — » oo. In order 
to work out the short- distance behavior, we demand that m{z) ^ nii when 
z is small (but still larger than microscopic scales). This is motivated by and 
consistent with the field-theoretic short-distance expansion [21,2], where, on 
a more fundamental level, field operators near a boundary are represented in 
terms of boundary operators multiplied by c-number functions. 

Now, the dependence of mi on hi is given by mi ~ /i^ " , and with the 
scahng relation = ^i/Pi we obtain that the scahng function MgpiO in 
(2.5) should behave as ~ C^i ^'^ in the limit C — 0. Inserting this in (2.5) leads 
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to the short-distance behavior 

m{z) - /i^''" ^(/^r-/?)/'^ for z-^lsp, (2.6) 

in agreement with Refs. [12,13]. In other words: Near the special transition 
a small hi gives rise to a macroscopic regime near the surface of depth Igp 
on which the behavior of the OP is governed by an exponent different from 
the one describing the long distance behavior for z ^ l^p. For the three- 
dimensional Ising model the short-distance exponent (/3^^ — /5)/z/ is negative, 
its (MC) literature value being —0.15 [17] (compared with —P/v ~ —0.52 [16] 
that governs the decay for z ^ Isp) ■ 

Next let us consider the case of large c, near the ordinary transition, the 
situation which is more natural if one thinks of possible applications for ex- 
periments and which was discussed in Ref. [18]. Since c is a so-called dangerous 
irrelevant variable it must not be simply set to oo from the start [2]. A careful 
analysis reveals that close to the ordinary transition the linear scaling field 
is given by hi = /ii/c^ with y = (A^^ — A^''^)/^. In MF theory the value of 
the exponent y = 1. Further, it is discussed in detail in Ref. [2] that in the 
framework of the e expansion one does not capture the deviation from the MF 
value in this exponent, while e.g. the 2;-dependence of expectation values is 
reproduced correctly. 

Analogous to (2.1), also near the ordinary transition the magnetization is a 
homogenous function of the linear scaling fields: 

m{z, T, hi) ~ r'^/^ m{z h-\ rb^/", hi 6^?^>) . (2.7) 

Removing b by setting it ~ 2; and setting t — again, one obtains the scaling 
form 

m{z, hi) ~ Z-^/" Mord(z/lord) , (2-8) 

where 

lard^K ' ' (2.9) 

is in this case the only length scale determined by the surface fields. 

In order to analyze the short- distance behavior of the magnetization, we de- 
mand again that m(2; ^ 0) ~ mi. Since for large c the surface is paramagnetic 
and responds linearly to a small external field, we have now mi ~ hi [14]. The 
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immediate consequence of the simple linear response for the scahng function 
AiordiC) occurring in (2.8) is that it has to behave as ~ (^^i"^ ^'^ for ( ^ 0. In- 
serting this in (2.7), we obtain that for z <^ lord the magnetization is described 

by 

m{z) ~ hi z'^ with k = ^ ^ ^ (2.10) 

Using that Ai/z/ = {d — r]\\)/2 and /3/iy = {d — 2 + 7])/2 together with the 
scahng relation ri± — {ri-\- ri\\)/2, one obtains k—1 — -qf^. 

In the MF approximation the result for k is zero. However, a positive value 
and, consequently, a non-monotonic profile are obtained when fluctuations are 
taken into account below the upper critical dimension d* = 4. The numeric 
value of K is 0.21, where we have used the Monte Carlo hterature value 0.8 for 
13^'^ /v [22] together with the scahng relation Ai + /3i = (d - 

2.2 Crossover between Special and Ordinary Transition 

In the previous section we discussed the behavior of the OP in situations 
where c — c* and c — > oo and where /ii 7^ 0. In both cases we found a 
near-surface behavior governed by the exponents of the special or ordinary 
transition, respectively, and a crossover to the power law z'^/" characteristic 
for the normal transition. 

However, the scahng form (2.2) should also describe the behavior of m{z) 
in the general case, in the presence of two length scales, and, in particular, 
it should cover the crossover between "special" and "ordinary" near-surface 
behavior. Let us assume that hi is fixed and vary c between c* and 00 such 
that the length scale Ic ranges from 00 and 0. First, in the limit Ic ^ 00 we 
were lead from (2.4) to (2.5). Due to the fact that the hnear scaling field at 
the ordinary transition is a combination of the surface parameters hi and c, 
the second limit, z/lc 00, is more subtle. In this limit the scaling function 
M-crit in (2.4) should behave as 

Mcrit{z/lsp, z/lc ^ 00) Mord({z/lspT ■ {z/lcf-'') , (2.11) 

where a = Af'/A"^'^ such that the product in the argument of Aiord on the 
right-hand side of (2.11) is just z/lgrd- 

Although experimentally a system where both scales lord and Ic are macro- 
scopic is probably not realizable, in the field-theoretic model to be analyzed 
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Z/Isp 

Fig. 1. Qualitative scenario for the crossover between special and ordinary transition 
for the order parameter. In our representation hi is fixed and c varies from c = c* 
to a value c S> c*. 

below there is no principle limitation. So at least from the theoretical point of 
view it is an interesting question to ask what shape the OP profiles have in the 
crossover regime between special and ordinary transition. In the following we 
present a qualitative scenario that, to our present knowledge, describes this 
crossover. 

Consider a system with fixed hi and variable Ic- As long as > 0, the asymp- 
totic behavior for small distances, z ^ Z^) is described by (2.6) where the 
amplitude of the short-distance power law does not depend on c. On the other 
hand, for distances large compared to all scales induced by hi and c, the OP 
profile should behave as ~ z~^l^ ^ again with an amplitude independent of c 
(and hi). More precisely, the crossover to the "normal" behavior takes place 
at a distance lint that varies between Ig^ for c — (f and lord for large c. 

The qualitative shape of the crossover profiles is depicted in Fig. 1. As long 
as > linti the OP is a monotonic function of with a single crossover at 
z = lint — hp- For If. < lint, howcvcr, the profile becomes non-monotonic and 
exhibits several crossovers. For z > Ic it first drops and thereafter increases 
again to approach the long-distance power law. For Ic ^ eventually the 
increase for z < Unt — lord must approach the power law (2.10). Of course, 
in this limit it would be appropriate to rescale the z-axis or readjust hi such 
that Icyrd remains finite. 
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2.3 The Correlation Function between Ordinary and Normal Transition 

As for the OP in the previous sections, a similar scahng analysis combined 
with arguments taken from short-distance expansions can also be carried out 
for the correlation function. Here we restrict ourself to a qualitative discus- 
sion of the two-point correlation function in directions parallel to the surface, 
called parallel correlation function in the following, for a system with large c. 
This function is useful to illustrate the general features one has to expect for 
correlation functions in the crossover regimes and later on it is needed for an 
heuristic argument that explains the near-surface growth expressed in (2.10). 
Further, we restrict the discussion again to the case r = 0. 



Fig. 2. Regimes of universal behavior (gray areas) of the correlation function in direc- 
tions parallel to the surface. The respective power laws governing the p-dependence 
of G{p, z) are shown. Crossover regimes are represented as white areas. 

The different regimes which the parallel correlation function G(p, -2, z) (where 
p denotes the distance between two point in a plane parallel to the surface) 
goes through are depicted in Fig. 2. First of all, for p < 2; the function G 
asymptotically behaves as in the bulk, G is governed by the bulk exponent 
Tj (lower light gray area). In the opposite case, p > z, the result depends on 
whether the distance from the surface is smaller or larger than the length 
scale lord introduced in (2.9). In the case z > lord for p ~ 2; the crossover from 
bulk to surface behavior takes place, where G is governed by the anomalous 
dimension of the normal transition (upper light gray area). The situation for 
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^ < lord is different. Again at p ~ 2; a crossover to surface behavior takes place, 
to a regime z < p < lord where the behavior is governed by exponents of the 
ordinary transition (dark gray area) . Eventually, for p ~ lord the crossover to 
the regime with "normal" behavior takes place. In Fig. 2 the crossover regions 
are represented as white areas. 

As we shall see in Sec. 3.3, this crossover behavior leads to the result that 
for instance in scattering with finite penetration depth of the radiation and 
lord greater than this depth, the structure function is effectively governed by 
exponents of the ordinary transition. 

2.4 Modifications atT ^ 

The scaling analysis presented above can be straightforwardly extended to the 
case r > 0. In d > 2, we may assume that the behavior near the surface for 
2; << ^ is unchanged compared to (2.6) and (2.10). The behavior farther away 
from the surface depends on the ratio l/^, where I stands for Isp or lord- In the 
case of / > ^ a crossover to an exponential decay will take place for -2 — ^ and 
the regime of nonlinear decay ~ z^'^^^ does not occur. For I < ^ & crossover 
to the power-law decay ~ z'^/" takes place for z ~ Z and finally at 2; ~ ^ the 
exponential behavior sets in. 

Similar results hold for the correlation functions. As long as for example the 
distance p in the parallel correlation function is smaller than the discussion 
of Sec. 2.3 holds. For distances larger than ^ correlations decay exponentially. 

An interesting phenomenon can be observed in the case ,^ < /. As said above, 
m{z) then never reaches the regime with power-law decay, but crosses over 
from the near-surface regime directly to the exponential decay. Let us restrict 
the further discussion to the case of the ordinary transition. Since the region 
where m{z) grows then only extends up to the distance ^, the maximum of 
the OP profile has roughly the value rrimax — hi where the exponent k was 
defined in (2.10). Now, the amplitude of the exponential decay should behave 
as ~ rrimax such that for 2; ^ ^ we have 

m(^)~hirexp(-^/0. (2.12) 

In other words, in the case ^ < lord the exponent k not only governs the behav- 
ior of m{z) near the surface, but also leaves its fingerprint much farther away 
of an universal dependence of the amplitude of the exponential decay. Nothing 
comparable occurs when ^ = 00 (compare Sec. 2.1 above), where all profiles 
approach the same curve m{z) ~ Az~^/'^ for z^lord, with an amplitude A 
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independent of hi. A similar phenomenon, termed "long-time memory" of the 
initial condition, does also occur in critical dynamics for T > Tj! [24]. The 
analogies between surface critical phenomena and critical dynamics will be 
discussed in more detail in the next section. 



2.5 Relation to Critical Dynamics and Heuristic Arguments 

The spatial variation of the magnetization discussed so far, especially the in- 
creasing profile at the ordinary transition, strongly resembles the time depen- 
dence of the OP in relaxational processes at the critical point. If a system with 
nonconserved OP (model A) is quenched for example from a high-temperature 
initial state to the critical point, with a small initial magnetization m^'\ the 
OP behaves as m ~ m^'^ t^ [25] , where the short-time exponent 9 is governed 
by the difference between the scaling dimensions of initial and equilibrium 
magnetization divided by the dynamic (equilibrium) exponent [26]. Like the 
exponent k in (2.10), the exponent 6 vanishes in MF theory, but becomes 
positive below d*. For example, its value in the three-dimensional Ising model 
with Glauber dynamics is 61 = 0.108 [27]. 

The high-tcmpcrature initial state of the relaxational process is to some ex- 
tent analogous to the surface that strongly disfavors the order and that (for 
hi — 0) belongs to the universality class of the ordinary transition. Further 
expanding this analogy, heating a system from a low-tempcrature (ordered) 
initial state to the critical point would be similar to the situation at the ex- 
traordinary transition. Eventually, analogous to the special transition would 
be a "relaxation process" that starts from an equilibrium state at Tj^. 

Motivated by the analogy between surface critical phenomena and critical 
dynamics, an heuristic argument for the short-distance growth of the OP at the 
ordinary transition behavior can be brought forward. In relaxation processes, 
after the system has passed through a non-universal regime immediately after 
the quench, the growth of correlation length ^(t) is described by a power 
law. Analogously one may argue that in the static case in the surface-near 
regime there exists an effective parallel correlation length ^11(2;) that grows 
as ^ z when the distance z grows away from the surface. As discussed in 
Sec. 2.3, for lateral distance p z the parallel correlation function crosses 
over from "bulk" to "surface decay", the latter being much faster than the 
former (compare Sec. 3.3). 

For the sake of the argument, consider the semi-infinite Ising system at bulk 
criticality where initially hi — with the tendency to order in the surface 
reduced compared with the bulk. If then a small hi is imposed on the sur- 
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face, a number of spins will change their orientation to generate a small mi. 
Now, the magnetization induced by in a distance z from the surface is deter- 
mined by several factors. First, it is proportional to nii and, in turn, for small 
mi, it is proportional to hi. Second, it should be proportional to the perpen- 
dicular correlation function C{0,z) ~ which says how much of a 
given surface configuration "survives" in a distance z. Eventually, it should 
be proportional to the correlated area in the distance z that is influenced by 
a single surface spin. According to the discussion above, the latter grows as 
^11 ~^ ~ z^~^. Taking these factors together we obtain 

m{z) ~ hi C(0, z) ~ hi z'-"!'"" , (2.13) 
in agreement with the short-distance law (2.10). 

The above argument holds as long as z and mi are small enough. If z is larger 
than the average distance between the spins flipped, the correlated areas start 
to overlap and the above argument breaks down. This is where the crossover 
to the normal transition sets in. 

From the fleld-theoretical point of view in both cases (surface critical phe- 
nomena and critical dynamics) the modified power laws are due to additional 
short-distance or short-time singularities near spatial or temporal boundaries. 
In both cases the singular behavior near the surface can be described by means 
of short-distance expansions or (in simplified form) by means of a phenomeno- 
logical scaling analysis as presented above. It is interesting to point out, how- 
ever, that in the presence of both spatial and temporal boundaries no addi- 
tional singularities (other than the ones discussed above) occur [28], and, as 
a consequence, the short-time behavior in critical dynamics near surfaces is 
governed by exponents that can be expressed in terms of known surface and 
dynamic exponents [28]. 



3 Mean-field (Zero-Loop) Approximation 

In this section we define the model and describe the results of the MF theory 
for general c and hi. They are the building blocks for the (one- loop) pertur- 
bative calculation of Sec. 4. Especially the zero-loop propagator for general c 
and hi is not contained in explicit form in the literature so far. Further in 
Sec. 3.3 we discuss the behavior of the real-space correlation function and the 
respective structure function that is measured for instance in X-ray scattering 
experiments. The one-loop calculation for this quantity for general c and hi is 
not attempted in the present work. So the MF results presented here are not 
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only the best presently available in the framework of continuum field theory, 
but also provide, as we think, a qualitative description of the results to be 
expected from scattering experiments. 



3.1 Model 



We consider the semi-infinite scalar 0^ with the Hamiltonian[2] 

1 



n 



V 



+ 



dV 



-ccf) +hi(p 



(3.1) 



where V and dV stand for the volume of the semi-infinite system and its 
(planar) surface, respectively. 

The well-known solution to the OP profile, minimizing the above Hamiltonian, 
takes the form[29,14,12] 

m^^\z) = a-\z + z+)-\ (3.2) 



where the length scale z~^ is given by 

, c + Vc^ + 4/iiq; 



(3.3) 

2ahi ^ ' 



and 



The solution (3.2) satisfies the general (mixed) boundary condition 

dm 

— — = cm — hi . (3.5) 
az 



The MF solution to the OP profile has a fairly simple form. Consistent with the 
discussion in Sees. 1 and 2, it is a monotonic function of z. Further, although 
c and hi in principle give rise to two independent length scales, in the OP 
profile they show up in a single length scale z'^ only. As a consequence, the 
shape of the profile does not change qualitatively when one moves from one 
fixed point to the other. 
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3.2 Free Propagator 



The two-point correlation function for general c and hi that will be described 
next solves the equation 



-V^ + -um'^{z) 



G'(x,x') = 5(x-x') , 



(3.6) 



subject to the boundary condition (at 2; = 0) 

d 



dz 



G'(x,x')|.=o = 0, 



(3.7) 



where m'^'^^z) is the profile of (3.2). To calculate G we Fourier transform with 
respect to the parallel, translationally invariant spatial directions. The result 
that can be obtained by standard methods reads 

G(p; z, z') = e{z - z')W{p, z)U{p, z') + e{z' - z)Wip, z')U{p, z) , (3.8) 



where p = |p|, 9{z) is the Heaviside step function, and the function W and U 
are given by 



W{p,z) 



2 , 3p 3 
P H : r + 



z + z+ {z + z+y 



-pz 



(3.9) 



and 



U{p,z) = — {B{p) 



2p5 



2 , 3p 
p + 



+ 



z + z+ (^ + z+Y 
2 3P 



-pz 



p 



+ 



with 



B{p)^ 



Z + Z+ (2; + z+Y 

-k'^{q + l) + 3kq-3q 
+ k^{q + l) + 3kq + 3q ' 



(3.10) 



(3.11) 



Further, we introduced a dimensionless wave vector 

k — pz'^ 



(3.12) 
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and, for the sake of brevity, we defined 

q = 2 + cz^ 



(3.13) 



The length scale z'^ was already introduced in (3.3). 

In contrast to the OP profile, the MF correlation function now depends on 
two length scales. In addition to the dependence on 2;"'", it contains via the 
quantity q defined in (3.13) an explicit c-dependence. The latter defines the 
second length scale l^. = c~^. 

For the correlation function G(p; 2;) and the for the propagator at p = 0, 
one obtains from (3.8) 



2 , 3p 3 

V ^ — \ — 7 + 



z + 2;+ (2; + z-^Y 



-2pz 



9 



{z + z+)^ {z + z+y 



(3.14) 



and 



G{0;z,z') = ^i [z + z^ 



X 



5 — q 
Q 



(z+)^ (z' + z+y^ + (z' + z+f 



9{z - z') ^{z^z') \ , (3.15) 



respectively. The results (3.14) and (3.15) are needed in Sec. 4 to calculate the 
one-loop contributions to the OP profile. 



3.3 Structure Function for Critical Diffuse Scattering 

A quantity directly relevant for experiments is the structure function 

00 CO 

S{p, K) = Aff G(p, z, z') e^('^^-«*^') dzdz' , (3.16) 







where A is a constant - for a cubic lattice it is just the illuminated area divided 
by a power of the lattice constant - and k — Kr-\-i Kiis a, complex wavenumber 
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taking into account that the radiation penetrates the material roughly up to 
the depth [4]. 

Analytic results for S on the basis of (3.8) can be obtained at or near the 
fixed points only, where the correlation function turns out to be much simpler 
than the general expression (3.8). Here we present the analytic results for 
the ordinary and the extraordinary (or normal) transition and evaluate the 
transformation (3.16) numerically in the crossover regime, i.e. for large c and 
and arbitrary hi/ c. The latter is the generic situation met in systems where the 
spin-spin interaction near the surface is not enhanced and where the surface 
spins arc coupled to an external (non-critical) medium that favors one of the 
spin directions. The results presented below were obtained by setting g = oo in 
(3.8) (thus ignoring corrections of order c~^). Then the function B{p) becomes 

lim B(p) = -— — , 3.17 

and the length scale z'^ from (3.3) is to be identified with the length scale lord 
introduced in Sec. 2.1. Again in leading order in c^^ one obtains 

hrd = Z'^ — —— , (3.18) 

consistent with the scaling analysis in Sec. 2.1 and with the fact that y = 1 in 
MF theory. 

By general arguments (scaling analysis and short-distance expansion) it was 
shown by Dietrich and Wagner [4] that in the limit p — > the structure function 
takes the form 

S{p,n)^so{n)-si{n)f\\-^ . (3.19) 

Explicit results were obtained in Ref. [4] at the ordinary transition up to one- 
loop order. 

Starting from (3.8) our result for the ordinary transition {t — hi — c — oo) 
is 



16 



in consistency with Ref. [4]. For the normal transition {t — 0, c — hi — oo) 
we obtain 



S{p, k) = 



20 k? 



2Fi(l,3,5,i«rV/2) 



12 

75 ' 



(3.21) 



where 2F1 is the hypergeometric function [30] and K denotes the real part 
of the term in square brackets. The MF values of r/y are 2 for the ordinary 
transition and 6 for the normal transition, respectively, and so the results 
(3.20) and (3.21) are consistent with (3.19). 

As the next step the transformation (3.16) is carried out numerically in the 
crossover regime between ordinary and normal transition. The result can cer- 
tainly not be compared with experiments on a quantitative level - the real 
exponents differ significantly from their MF values - , but it should give the 
qualitative picture that one has to expect from a (Gedanken)experiment car- 
ried out in a system at bulk criticality with adjustable scaling field hi/c and 
otherwise fixed parameters. 

0.25 

S lord 
0.2 

0.15 

0.1 

0.05 





0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

p lord 

Fig. 3. The MF structure function for critical diffuse scattering defined in (3.16) in 
the crosover regime between ordinary and normal transition for different values of 
lord Ki (solid lines) . The upper dashed line shows the asymptotic {p —>■ 0) result at 
the ordinary transition given in (3.20). The lower dashed line is the analogous curve 
at the normal transition of (3.21). 
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The results for S for the particular choice depicted in Fig. 3 

(solid lines). For / = lord ' i^i 3> 1, the p ^ limit of S has effectively the 
cusp shape typical for ordinary transition. (In contrast to the situation in the 
bulk, S never diverges[4], because in the surface- near regime (which at Tj! 
extends to infinity) the long-distance decay of correlations is much faster than 
in the bulk.) For smaller values of it becomes visible that for small momenta 
plord ^ 1 the behavior is actually described by the anomalous dimension of 
the normal transition, leading to a curve without cusp departing with slope 
zero from p = 0. For comparison we have also plotted the asymptotic results 
for p — > and lord — and oo (dashed lines) . 



4 One-Loop Calculation of the Order-Pcirameter Profile 

In this section we calculate the OP profile in the framework of the e-expansion 
(where e = 4 — d) near both the ordinary and special fixed point, the latter 
mainly with regard to the crossover between special and ordinary transition. 
For the sake of simplicity, all calculations are restricted to r = 0. In the 
following we first give the results for the one- loop contribution to m{z). They 
are regularized dimensionally in d — A — e dimensions and renormalized by 
minimal subtraction. After focussing attention on the short-distance behavior 
of the OP, we calculate the full scaling functions up to first order in e near 
both the ordinary and the special transition. Finally we present results on the 
crossover between special and ordinary near-surface behavior. 

In order to remove the UV divergences appearing in the one-loop term, the 
parameters u, c, and hi of Sec. 3 have to be considered as bare parameters 
labelled with an index '0' below. In the standard procedure, dimensional reg- 
ularization combined with minimal subtraction [2], the bare parameters are 
replaced by renormalized ones to obtain a UV finite expression for m{z). The 
field (f) which in general also has to be renormalized is unchanged up to order 
e, such that we do not have to worry about wavcfunction renormalization in 
the present work. Further, in dimensional regularization the fixed point value 
c* of the surface enhancement remains at zero such that c — c. 

4-1 Perturbation Expansion 

Expanding in terms of the coupling constant uq, m{z) can be written in the 
form 

m{z) = m^'^) (z) + m« (z) + 0{ul) , (4.1) 
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where the MF solution m^^^ was given in (3.2) and the first-order term m^^^ is 

oo 

= -y / dz'G(0,z,z') J G(p, z' , z') m^''\z') (4.2) 



with Jp = (27r)-''+i / d'^-V- Using the results (3.2), (3.14), and (3.15), the z' 
integration in (4.2) can be carried out. This yields the result 



(1) 



^+r((3-e)/2) 



a-^X J dkk-'I 



(4.3) 



where Sd — (2-\/7r) and the integrand / is given by 



I ^-BFe-^'''/'^ + -X 

2 



k-^ + —XB 
2q 



+^ [B X{q + 9) + 3X{q - 1) - 3g] k'' 
fX^ 1 ^ /5 3\ , , 



(4.4) 



where 5 and k were defined in (3.13) and (3.12), respectively. Further, we 
introduced 



X 



z + Z' 



+ 



(4.5) 



with z~^ defined in (3.3) and 



F = -k-' + -Xk-^ + -X^k-^ . 
4 4 2 



(4.6) 



4-2 Renormalization of UV Divergences and Short-Distance Singularities 



The perturbative result (4.3) contains UV singularities which show up as poles 
oc e"^ in dimensional regularization. These poles are removed by minimal 
subtraction[2,31], by expressing the bare parameters occurring in m^^^ in terms 
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of renormalized ones. The relation between bare and renormalized coupling is 
given by 



UoSd = 1^ u 



1 + ^ + Oiu") 



(4.7) 



where ^ is an arbitrary momentum set to 1 afterwards. Further, near the 
special transition, for finite c and /ii, the renormalizations are 



Co = 



1 + ^ + Oiu') 



(4.8) 



and 



u 



1 + ^ + 0{u') 



(4.9) 



At the ordinary transition, for c — > oo with finite hi — hi/c, we have 



hi ,0 = //'hi 



u 



1 + ^ + 0{u') 



(4.10) 



Up to first order in e the renormalization of hi and hi are the same, but taking 
into account higher orders in e they are different [2]. In both cases the poles 
are removed and one is left with an expression for ■m{z) that is regular in e 
and u. 

In our first-order perturbative calculation, modifications of the short- distance 
behavior of the OP profile show up in the form of logarithmic (short-distance) 
singularities of m^^\ The source of these logarithms in our result (4.3) is the 
first term on the right-hand side of (4.4). Focussing the attention for the 
moment on the short-distance behavior only, (4.3) can be written as 



^(^) ~ "^"^ — Jsing + terms regular for ^ — >• , 



(4.11) 



where Jging is given by 



OO 

J sing ~ J dk B k ^ e ^^/^ 



(4.12) 
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For a closer inspection of the z ^ limit of this integral, it is useful to rewrite 
the function B (defined in (3.11)) in the equivalent forms 

B^l- 2bsp = -1 + 2kbord (4.13) 



with 



b = + (414) 

k^ + {q + l)k^ + 3qk + 3q ^ ^ 



and 

k^ + 3q 



-'ord 



k^ + {q + 1)F + 3qk + 3q 



(4.15) 



While bgp yields regular terms for general c and hi, in the second form bord 
leads to finite contributions after taking the limit c — > oo: 

3 

lim bord = &00 = ,0 , o7 I o ■ ('^■16) 



With (4.14) and (4.15), Jging can be written as 

Jsing^ J dk{2kbord - l)k-'-'e-^''/'^+ J dk{l - 2bsp)k''-'e~^'''/'^, (4.17) 



where we assumed that cz^ > 1. Alternatively, after straightforward manipu- 
lations, it can be expressed in the form 

oo 

Jsing = j dkk-^-' (2e-^'''/^- - e"^*^^/^^) + terms regular for ^ ^ , (4.18) 
1 



where lc = c ^ is the length scale induced by the surface enhancement. 

To discuss the crossover between special and ordinary near-surface behavior, 
we assume c to be large and /c ^ — ^ord- From (4.18) it becomes then 
apparent that for z <^lc 

J sing — log-z + terms regular for 2; — > . (4.19) 
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For /c -C ^ -C z^^ on the other hand, the first term in the integral in (4.18) is 
neghgible, and we obtain 

Jsing — — log z + terms regular for 2; — > , (4.20) 

which also holds for c — > 00 and z <^ z'^ , i.e., it describes the short-distance 
singularity at the ordinary transition. 

If eventually the coupling constant u is set to its fixed-point value u* = e/3 
and the logarithms are exponentiated, we find 

[ 1 - I log z + 0{e^) ~ z-'/^ for z^h 
m{z) ~ J . (4.21) 

[ 1 + I log ^ + 0(e2) ~ for Ic^Z^ lord 

The exponents e/6 and — e/6 arc to be identified with the first-order values 
for 1 — rjf^ and (/^^^ — j3)/v^ respectively [2]. In general, for finite cz'^ both 
terms in (4.17) give non-negligible contributions, each on the respective length 
scale, i.e. the first term only for z <^ the second one for z <^ lord- This is 
consistent with our qualitative discussion in Sec. 2.2. 



4-3 Scaling Function at the Ordinary Transition. 



The procedures for obtaining the scaling function at the ordinary and the 
special transition are similar. We describe the first case in more detail. For c 
00, after the subtraction of pole terms with (4.7) and (4.10), the perturbative 
results (4.1) and (4.3) take the form 



m{z) — mP{z; u, hi) — 



2usdX 



z^a 



00 

J dkC'"''^ + ^{2X -3) log z^ 
.0 



, (4.22) 



where 



C-<^ = F(1 - 6oo)e-'"^/"^ + ^X{b^ - X)k-' + ^X{b^ - l)k 



-2 



+ j(2^ - 3) 



(4.23) 



is obtained from I (defined in (4.4)) by setting q — c — 00. In dimensional 
regularization terms of order 0(A;°) and 0{k^) present in (4.4) give vanishing 
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contributions to the integral and are omitted here. As the next step we utihze 
the scahng form (2.8), derived on a more rigorous level with the help of the 
renormalization group [2]. The scaling function is given by 



MordiC) = m(z = 1, = u*, hi(C) = hiZ^?"''/-) , 



(4.24) 



where the scaled distance ( is 



C = 7^ with lord = (a* hi) ''/^i'''' , 

I'ord 



(4.25) 



and a* is defined by (3.4) with u — u*. 

The result (4.22) gives the scaling function up to terms O(e^) according to the 
equation 

MordiO = m^''\l,u\ hi^^rV^) hiz^?"/'^) + 0{e') (4.26) 



The MF result (3.2), after substituting 2; — > 1 and 

'h^)-'z-^r^/^ = (1 + I logC) + 0{e') (4.27) 



z (a 



reads 



m 



(0) 



Together with the one- loop contribution (4.22), we obtain for the scaling func- 
tion 



-1 c 



1 + C 



1-6 l^iogc+ydkc-''' 



+ 0{e^) . (4.29) 



One can easily check that for — >• the integral /q°° dk C"''"'^ behaves as ~ 
— |logC+ regular terms, and, as a consequence, the short-distance limit of 
the scaling function is given by 



MordiC) ~ C 



l-glOgC 



+ 0{^) . 



(4.30) 
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Taking into account that = 1 - e/2 + 0{e^) such that z"^/^ = ^"^(1 + 
e/2 log 2; + (9(e^)), this is consistent with (4.21) and, hence, with the result 
(2.10) of our scaling analysis. Furthermore, it is straightforward to show that 
■M-ord approaches a constant for ^ — > 00. 

In order to exponentiate the scaling function in the short-distance regime, we 
effectively replace the logarithm of (4.30) by the power law "^^'^i the exponent 
1 — e/3 being the one-loop value of A.'^'^ /u. This is achieved by adding 



yard 6 



l-e/3 



^0 



glOgC 



(4.31) 



to (4.29), a function of 0{e^). 

The result for the OP as a function of i.e. the function 

Mord = C'^'^^Sord + a'-Mord) , 



(4.32) 



is depicted in double-logarithmic representation for three different values of e 
in Fig. 4. The integral in (4.29) had to be calculated numerically. All profiles 
show the characteristic crossover between short-distance growth and long- 
distance decay. However, extrapolating e to one obviously leaves us with a 
profile with somewhat unexpected features. It has a point of inflection were 




0.001 



100 1000 



Fig. 4. Scaling function Mord for the crossover between ordinary and normal 
transtion defined in (4.32) for e = 1 (solid), 0.4 (dashed), and 0.1 (dotted line). 
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one actually expects a smoother scaling function. On the other hand, for the 
smaller values of e this problem does not exist. So wc think that extrapolating 
e to 1 in (4.29) probably just at the limit, where the first-order result ceases 
to provide a good approximation. 



4.4 Scaling Function at the Special Transition and Crossover to the Ordinary 
Transition 



For arbitrary c wc proceed as before, now with the rcnormalization of the 
surface parameters given by (4.8) and (4.9). The perturbative result for the 
profile (4.1) is [13] 



m{z) — m^^\z; u, c, hi) 



2uSdX 



,(4.33) 



where we extracted the regular contribution 



C = {2b,p -1)F e-2*^^/^+ + -X 



+ 



{q + 9)X 
2q 

1 fx 3 



2\q 2) k{k + l) 
from / of (4.4). 



, q + l\X 



(4.34) 



At the special transition, alternatively to (2.2), the order parameter can be 
expressed in the scaling form 



m{z]u, c,^i) = ^-^/"M(C,7) 



(4.35) 



with 



M{C, 7) = m{z = l-u = u\ c(C, 7) = cz^'\ HO = /^i^^"/^) 
= m(°) {l;u*, cz''/\ hiz""'"'^) + m^^\l- u*, cz^/\ hiz""'" + 0(e')(4.36) 

where 



C = 1— , J = chi , and Igp = {a* hi] 

''sp 



(4.37) 
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Expressing all the parameters in terms of scaling variables ( and 7 and keeping 
terms to first order in e in the MF contribution and to zeroth order in e in the 
one-loop term, we obtain the result 



1 Co 



log Co 



+ 



|logm+2/d^C(^) 



+ O(e') , (4.38) 



where 



Co= 2 (7 + \/7' + 4 



(4.39) 



and the integrand C is given by (4.34) with z — 1, — Co/C? Q — 2 + Co7- 

It is straightforward to verify that the short-distance behavior of the OP 

described by (4.38) for a system with 7 ^ 1 is consistent with the discussion 
in Sec. 4.2, especially with (4.21). The short-distance logarithm in (4.38) is 
caused by the first term on the right-hand side of (4.34). One obtains 



A<(C,7)~f 



(4.40) 



for C/^o ^ Q ^ (equivalent to z <^ Ic), which, together with the prefactor 
~ C~^^^ (see (4.35)), is consistent with the upper line in (4.21) and with 
Refs. [12,13]. And one finds 



A^(C,7) 



Co 



(4.41) 



for g ^ <^ C/Co ^ 1 (equivalent to Ic <^ z <^ lord), in agreement with the 
second case in (4.21). 

The crossover between the two asymptotic cases occurs at the length scale Igp, 
which is the intermediate scale between Ic and lord- The ratio between all the 
scales is set by the single parameter 7 of (4.37), and from the definitions of 
the length scales (2.3), (2.9) and 7 we find 

lord isp 

i_^^WAr ^ = 7^/^ (4.42) 
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At the special transition, 7 = 0, the scahng function is given by Aigp — 
-^(CjT = 0)1 with A4 given by (4.38). Like at the ordinary transition we 
exponentiate the short-distance logarithm by adding the C(e^) function 



^i-V3_^(l_ZelogC 



(4.43) 



where the exponent 1 — 2e/3 is the first-order value for (31^ /v. The OP scahng 
function 



Msp = C"'+'/' {Qsp + a*Msp) 



(4.44) 



is depicted in Fig. 5 for three values of e. All profiles correctly show the asymp- 
totic power laws for short and long distances, respectively. Like the result for 
the ordinary transition, however, the curve for e = 1 exhibits unexpected fea- 
tures. This becomes clearer visible, when one plots the scahng function M.sp 
directly, without exponentiation, i.e. the result (4.38) with 7 = 1. The result 
is shown in Fig. 6 (upper curve). For z — 0, J^sp behaves as ~ logC, and 
for — cxD it approaches a constant. As seen from Fig. 6, it varies between 
the two asymptotic limits in a non-monotonic fashion, where one actually ex- 
pects a monotonic variation. The reason for this qualitative defect was already 
discussed at the end of Sec. 4.3. 




0.001 0.01 



10 100 
Z/Isp 



1000 



Fig. 5. Scaling function M.sp for the crossover between special and normal transition 
defined in (4.44) for e = 1 (solid), 0.4 (dashed), and 0.1 (dotted line). 
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0.001 0.01 0.1 1 10 100 1000 



Fig. 6. The (unexponeniated) scaling function M. of (4.38) for different values of 7 
and e = 1. 

What happens if 7 takes intermediate values, between the fixed points 7 = 
(special) and 7 = 00 (ordinary)? The results for ^^((^,7) of (4.38) for 7 = 0, 
1, 5, 10, and 15 and e — 1 are displayed in Fig. 6 (from top to bottom). 
For increasing 7 the profiles first become monotonic functions and then, for 
7 >1, a minimum occurs. This is the signature of the crossover between special 
and ordinary transition. Consistent with the qualitative discussion in Sec. 2.2, 
drops first and then increases again to approach a constant. However, for 
larger values 7 ^15 the scaling function M. becomes negative, signaling that 
the one-loop result calculated at the special transition becomes meaningless 
at intermediate values of C,. This time the problem is not cured by going to 
smaller e; for large 7 the scaling function M. becomes still negative. 

Eventually in Fig. 7 we have plotted the function 

M^C^^''''{gsp + a*M) , (4.45) 

where Qsp is defined in (4.43), for the same values of 7 as in Fig. 6. All curves 
have the expected power-law form for small and large arguments. For inter- 
mediate values up to 7 ~ 10, the behavior is reminiscent to the qualitative 
picture described in Sec. 2.2 and depicted in Fig. 1. For larger 7^15 the OP 
becomes negative for intermediate z/lgp- 
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0.01 0.1 1 10 100 1000 



Z/Isp 

Fig. 7. Scaling function M for the order parameter in the crossover regime between 
special and normal transtion defined in (4.45) for different values of 7 and e = 1. 

5 Summciry and Concluding Remcirks 

We studied a semi-infinite Ising-like system by means of continuum field the- 
ory. Particular attention was paid to situations where the surface enhance- 
ment c and the magnetic field hi are off their fixed-point values, in which 
case they may give rise to macroscopic length scales (in addition to ^) and 
anomalous behavior of thermodynamic quantities near the surface. We car- 
ried out a phenomenological scaling analysis and discussed the near-surface 
phenomena from various points of view, especially analogies with in critical 
dynamics were pointed out. Furthermore, we presented mean-field results for 
the structure function in the crossover regime between ordinary and normal 
transition, and we carried out a first-order calculation for the order-parameter 
in the framework of renormalization-group perturbation theory (e-expansion) . 

Let us conclude by highlighting the main results of this work. The structure 
functions displayed in Fig. 3 may, on a qualitative level, directly be compared 
to the outcome of scattering experiments, like the one carried out by Mailander 
et al. [5]. Our results show that in the case when the length lord is larger than 
the penetration depth of the radiation k„^^, the structure function is effectively 
governed by the exponent of the ordinary transition and shows the cusp for 
p — > predicted in Ref. [4]. Only when the product lord ■ i^i is smaller, the 
signature of the normal transition, a "flat" structure function, can be seen. 
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Further, the ordcr-paramctcr profiles for the crossover between ordinary and 
normal transition displayed in Fig. 4 supplement earlier studies which focussed 
on various other aspects of the near-surface behavior [18-20]. Although extrap- 
olating e to one in our results appeared to be rather courageous, we believe 
that with these profiles a quantitative comparison with experimental data 
should be feasible. The Monte Carlo data available at present [19], although 
obtained from relatively large systems with up to 256^ x 512 spins, are still 
distorted by large finite-size effects and can not yet compete with the results 
presented in this work, where really the situation in the semi-infinite system 
is described. 

Eventually, the crossover between special and ordinary transition was studied 
on the basis of the first-order perturbative approximation. In principle the 
results for general c and hi presented in Sec. 4.4 should cover this crossover as 
the limit c oo. The qualitative scenario was described in Sec. 2.2. In our first- 
order calculation we found that, while the asymptotics for small and large z 
remain correct (compare also (4.21), the one-loop contribution becomes larger 
than the zero- loop term for intermediate z, and the order-parameter profile 
becomes even negative for large c (see Figs. 6 and 7). This means, in other 
words, that the first-order approximation breaks down in this limit, at least for 
intermediate distances, where the crossover from z^^^i -P)/^ (characteristic for 
the special transition) to z^^"'' -P)/'^ (characteristic for the ordinary transition) 
should occur. The reason for the failure of the first-order approximation is 
connected to different UV divergences at the ordinary and special transition, 
which, in the framework of the e-expansion, are intimately related to the short- 
distance singularities. This problem could probably be solved by introducing a 
cutoff, which, from the physical point of view, would certainly be meaningful, 
and consider the ordinary transition as the case where the surface enhancement 
c is of the order of the cutoff (as opposed to being of the order of the finite 
momentum scale in dimensional regularization) . 
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